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Abstract 

We describe a new cascade mass reconstruction approach to allow 
reconstruction of SUSY particle masses in long cascades (five or more 
particles) at the LHC with integrated luminosity as low as a few fb -1 . 
This approach is based on a consecutive use of the endpoint method, 
an event filter and a combinatorial mass reconstruction method. The 
endpoint method gives a preliminary estimate of light sparticle masses. 
An event filter combining the maximum likelihood distributions for all 
events in the data sample allows suppression of backgrounds and gives 
a preliminary estimate of heavy sparticle masses. Finally, SUSY parti- 
cle masses are reconstructed by a search for a maximum of a combined 
likelihood function constructed for each possible combination of five 
events in the data sample. 

SUSY data sample sets for the SU3 model point containing 80k 
events each were generated, corresponding to an integrated luminos- 
ity of 4.2 fb -1 . These events were passed through the AcerDET de- 
tector simulator, which parametrized the response of a detector. To 
demonstrate the stability and precision of the approach five different 
80k event data sets were considered. Masses were reconstructed with 
a precision of about 10% for heavy sparticles and 10 — 20% for light 
sparticles. 
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I. Introduction 



Supersymmetry [1] is a very popular extension of the Standard Model. By 
introducing an equal number of fermion and boson degrees of freedom SUSY 
leads to a cancellation of quadratic divergences in the Higgs sector and pro- 
vides an attractive solution to the hierarchy problem [2]. If supersymmetry 
exists at an energy scale of 1 TeV, the production cross sections of SUSY 
particles can be significant and effects of SUSY particles should be observed 
at the LHC, which should reach an integrated luminosity of 300 fb~ x in about 
5 years [3]. In many supersymmetric models, in order to avoid undesirable 
weak scale proton decay, R-parity conservation is assumed. This leads to 
pair production of SUSY particles and therefore the lightest SUSY particle 
(LSP) is stable. As a result, SUSY events at an accelerator will give two 
decay chains, each containing one LSP and Standard Model particles in the 
final state. The reconstruction of a SUSY event is complicated because of 
escaping LSPs and many complex and competing decay modes. In this paper 
we will consider mass reconstruction of five SUSY particles in the cascade 
decay 

g ^bb 2 ^ -> W2&1&2 Xihkbib 2 (1) 

The gluino decay chain (1) is shown in Fig.(l) 




Figure 1: A gluino cascade decay chain. 

At present there are two different approaches to SUSY mass reconstruc- 
tion. The standard technique for analysis is "the endpoint method" which 
has been widely studied [3] - [10] for high integrated luminosity of about 
100 — 300 fb _1 at the LHC. This method looks for kinematic endpoints. 
SUSY particle masses are reconstructed by minimizing the \ 2 function of 
the difference between observed and theoretical endpoints. 

The second method of SUSY particle mass reconstruction is the so called 
"mass relation method" [11], [12] based on the mass relation equation (see 
Appendix A) which relates SUSY particle masses and measured momenta of 
detected particles. The "mass relation equation" is obtained as a solution of 
a system of four-momentum constraints for each vertex in the decay chain 
(1). For example for the gluino decay vertex it follows that m~ = (p^o + 
hi + h 2 + Ph + Vb 2 ) 2 which give the relation between four components of 
LSP four-momentum in terms of decaying particle mass, LSP mass and the 
four-momenta of detectable particles. Similar relations can be obtained for 



each vertex in the process (1). There are four vertices in the cascade (1) so 
one gets four kinematic equations that can be solved for four components of 
LSP four-momentum in terms of SUSY masses and momenta of detectable 
particles. By substituting these components into the on-shell mass condition 
for LSP m|o = p|o the "mass relation equation" including all SUSY masses 
(m) and detectable particle momenta (p) in the process (1) 

f(m,p) = (2) 

is found (an explicit form of this equation is given in Appendix A). This 
equation is valid for any SUSY event including the cascade (1). For one 
event this equation is underconstrained, it contains five unknown masses. 

Authors of [12] developed two approaches to the problem of mass search 
at high luminosity of about 300 fb" 1 . Both approaches assumed that three 
light SUSY particle masses (x° 2) l R) and xT) are known. In this case, the mass 
relation equation contains only two unknown parameters, the masses of the 
heaviest sparticles. In the first "event pair analysis" approach it was assumed 
that SUSY particles are the same from event to event. Gluino and bottom 
squark masses are reconstructed in this case by solving the system of two mass 
relation equations. The second approach is based on the search for the max- 
imum of the combined likelihood function for all events which is constructed 
by an individual likelihood function for each event taking into account the 
mass relation equation in two dimensional (m^m^) mass space. The search 
for the maximum of the combined likelihood function allows reconstruction 
of gluino and bottom squark masses at high integrated luminosity. 

We note that the reconstruction at low integrated luminosity (with a 
few fb _1 ) does not work for heavy sparticle masses for either "the endpoint 
method" and "the mass relation method" in their original forms. 

"The endpoint method" uses only events near the end point; therefore, 
high statistics is required. Also, in order to reconstruct the gluino mass a 
non-relativistic approximation is used for the LSP near the end point. This 
approximation is not justified for situations with In close in mass either to 
X% or Xi- Note that in [9] an approach was proposed to obtain the gluino 
mass without a non-relativistic approximation. 

"Pair analysis" of "the mass relation method" assumes that masses of 
SUSY particles are the same from event to event while in reality heavy par- 
ticle masses follow a Breit-Wigner distribution. Therefore, an attempt to 
extend this approach to five masses often leads to an inconsistent system of 
equations or wrong solutions. An attempt to extend the original combined 
likelihood function approach [12] based on a grid search in five dimensions 
leads to unreasonably intensive computing calculations. 

The goal of the present article is to develop an approach allowing recon- 
struction all of five masses of the cascade (1) at low integrated luminosity of 
a few fb -1 . This luminosity can be reached at the early stage of the LHC in 
comparison with a projected integrated luminosity of 300 fb -1 in five years. 



This approach is based on a consecutive use of the endpoint method, an 
event filter and a combinatorial mass reconstruction method. It was found 
that reliable results are obtained only with a good starting point for the final 
fit and by rejecting a large fraction of the background. The first two methods 
are used to define five SUSY particle mass ranges for the final combinatorial 
mass reconstruction method in which all five masses are simultaneously fit 
and allowed to vary from event to event. The endpoint method therefore is 
used to get a preliminary estimation of the light SUSY particle (x^ Ir, Xi) 
masses and corresponding errors. These masses are then used with the mass 
relation equation constraint to construct the maximum likelihood distribu- 
tion in the two heaviest-sparticle mass-plane for each individual event. An 
event filter combining the maximum likelihood distributions for all events in 
the data sample allows determination of the range of heavy (g, b) masses and 
significant suppression of the background. Finally SUSY particle masses are 
reconstructed by a search for a maximum of a combined likelihood function, 
which depends on all five sparticle masses, constructed for each possible com- 
bination of five events in the data sample. The approach is self-contained; 
no information on the nature of the data set (SU3, SPSla or other) is used. 

Note that a possible way to improve "the mass relation method" was 
proposed in [13] by considering all events at the same time and fitting all five 
masses simultaneously. 

II. Simulation 

We choose for this study the SU3 model point. This point has a significant 
production cross section for the chain (1): gluinos and squarks should be pro- 
duced abundantly at the LHC. The bulk point SU3 is the official benchmark 
point of the ATLAS collaboration and it is in agreement with the recent 
precision WMAP data [14]. This model point is described by the set of 
mSUGRA [15] parameters given in Table (1). 
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SU3 
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300 GeV 


-300 GeV 


6 


> 



Table 1: mSUGRA parameters for the SU3 point. 



Branching ratios for the gluino decay chain (1) at the SU3 point are 

- 16^ ~ h 2U% ~0 11^% Ir 100% „o ^ Q 46% 

g^b^xV^lR^xi =► 0.18% 
Assumed theoretical masses of SUSY particles in the cascade (1), the total 
branching ratio and a cross section generated by ISA JET 7.74 are given in 
Table (2). 
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a\pb] 


SU3 


720.16 


605.93 


642.00 


223.27 


154.63 


118.83 


0.64% 


19 



Table 2: Assumed theoretical masses of sparticles, branching ratio BR and 
production cross section a at the SU3 point. Masses are given in GeV. 



Monte Carlo simulations of SUSY production at model points were 
performed by the HERWIG 6.510 event generator [16]. The produced 
events were passed through the AcerDET detector simulation [17], which 
parametrized the response of a detector. Samples of 80k SUSY events were 
used. This approximately corresponds to 4.2 fb _1 of integrated luminosity 
because the SUSY production cross section is 19 pb at the SU3 point. Five 
different sets of 80k SUSY events were considered to demonstrate the stability 
and precision of the mass reconstruction approach. 

In order to isolate the chain (1) the following standard cuts were applied: 

• two isolated opposite-sign same-flavor (OSSF) leptons (not tau leptons) 
satisfying transverse momentum cuts pr(^ ± ) > 20 GeV and Pt{P) > 10 GeV 

• two b-tagged jets with p T > 50 GeV; 

• At least three jets, the hardest satisfying p T1 > 150 GeV, pr2 > 
100 GeV, p T3 > 50 GeV; 

• M eff > 600 GeV and E™ iss > 0.2M e// , where E™ iss is the missing 
transverse energy and M e ff is the scalar sum of the missing transverse energy 
and the transverse momenta of the four hardest jets; 

• lepton invariant mass 50 GeV < M u < 105 GeV. 

Note that at the first stage of the reconstruction procedure to estimate 
light masses (xfj, ^R? an d X?), the chain q~ L — > x\q — > IrI^ — > Xihhq wa s 
considered with the same cuts as mentioned above except the requirement of 
b-tagged jets and without a cut in lepton invariant mass. Lepton invariant 
mass Mu cuts were determined from analysis of // distribution for the above 
q~L chain which gives the upper edge with relatively high precision. Mu cuts 
reject about 1/4 of signal events but improve the signal to background ratio. 

It was shown in [7] that the QCD processes are cut down by the require- 
ment of two leptons and of considerable missing The processes involving 
Z and W are suppressed by the requirement of high hadronic activity together 
with high missing E T . The only Standard Model background surviving the 
hard cuts mentioned above is tt production, where both W's decay leptoni- 
cally into a bbll state. The tt background contributes about 20% to the total 
number of events with bbl + l~ in the final state after cuts. We found that this 
number can be reduced to about 2% by applying an event filter procedure 
described in chapter IV. This contribution is negligible in comparison with 
more than 100% contribution due to SUSY bacground and we do not include 
tt background in the foolowing analysis. 

Table (3) shows the number of signal events and SUSY background events 



for the SU3 model point after cuts were applied to the five sets of 80k SUSY 
events that corresponds to an integrated luminosity of 4.21b" 1 . A classi- 
fication of events as signal and SUSY background ones is based on truth 
information. The SUSY background to the process (1) is significant. It fol- 
lows from Table (3) that for the SU3 point the number of SUSY background 
events is a factor 2 greater than the number of signal events. 
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3.3 
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50 
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3.1 
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141 


55 


86 


2.6 


1-5 


731 


244 


487 


3.0 



Table 3: The number of signal and SUSY background events after cuts ap- 
plied to 80k SUSY events. Ratio = (Signal+Background) /Signal. Row 1-5 
shows the number of events for five sets combined. 

A better suppression of the SUSY background could be achieved from de- 
tail examination of background contributions. Three different types of SUSY 
background can be considered: a background containing r leptons decaying 
into electrons or muons, a background with related processes in which \2 de- 
caying into Ir appears in a wrong vertex or if no gluino is produced as shown 
in Fig. (2) and other background. We estimate their relatives contributions 
at the SU3 point as 50%, 25% and 25 %, respectively. 

A significant contribution to SUSY background comes from dominant 
modes of g, b, and X2 decays. For the SU3 point xX decay into the pair 
rfi is favored over the decay to Ir due to the relatively high value of tanf3: 
BR(x° — ► T ± fi =F — » r ± r =F Xi) equals 48.7%. The presence of r leptons abun- 
dantly produced in dominant decay modes of SUSY particles would be a clear 
indication of SUSY background, r-decays can create electrons or muons in 
the final state thus faking the cascade (1). The following decay sequences 
of gluino or bottom squark also imitates the cascade (1) after r-decays into 
leptons: 

g^—ttxt—? bxibW + — > bf{~ u T bl + v { — > for" ~x\i> T bV ~v x , 
b\ -> xlt -> fri7 T bW + -> T'xlv T bl + vi . 

The application of r tagging could provide an additional suppression of 
combinatorial SUSY background by approximately a factor of 2. Note that 
in this work we did not use r tagging for background suppression. 

For the second type of background (Figure 2) only one b-quark appears 
in a wrong vertex which effectively leads to smearing in momentum of the 
corresponding b-jet. As a result, for reasonably small smearing these pro- 
cesses, even though they are background, can give quite correct mass peak 
positions upon reconstruction. 



Figure 2: Examples of SUSY background processes: neutralinos appear in 
wrong vertices (left); gluino is missing (right). 

III. Preliminary estimate of light masses by 
the endpoint method 

The endpoint method has been widely used to determine masses of SUSY 
particles [5] in particular in application to a decay chain 

QL X2Q ~^ llllnearQ ~^ Xlhar^nearQ (3) 

which is a subprocess of the cascade (1) if one considers q~z instead of b. For 
leptons in the process (3) the following notations are used: l nea r stands for 
a lepton which is "nearer" to the quark while the lepton radiated by the 
slepton is called lf ar because it is "further" from the quark. 

In the process (3) there are three detectable particles: a quark and two 
leptons. Four invariant mass distributions can be formed from them: mn, 
m qii, m qi n ear! m qi fa r • The m n an d m qii distributions are observable, in partic- 
ular, the 11 edge, qll edge and qll threshold can be measured. It is impossible 
to observe m q i near and m q i far distributions because the l nea r and lf ar assign- 
ment is ambiguous, however these invariant masses can be combined to give 
ql low and ql high edges which are observable. 

Endpoint extraction 

To extract endpoints from these distributions a fit was performed. The 11 
distribution was fitted in a standard way by a triangle form convoluted with a 
Gaussian. Usually the four remaining distributions are fitted near the edges 
by straight lines, but for 80k data samples this approach does not work, due 
to low statistics. Instead predefined functions covering a wider range of the 
distributions (0-1 TeV) have been used to get the endpoints. The qll, ql low 
and ql high distributions were fitted by a parabola and a straight line in the 
range 0-1 TeV. The qll threshold distribution was fitted by an exponential 
function and a straight line in the range 0-1 TeV. 
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Figure 3: Invariant mass distributions for Set 1 at the SU3 point: 11 edge 
(left), qll edge (right) 



SUSY Hist Edge 3 



] 



Entries 


316 


Mean 


474.3 


RMS 


163.3 


Left 


265.1 


Mean 


381.8 


Decay 


143.1 


Max 


26.4 



SUSY Hist Edge 4 



] 




Entries 


599 


Mean 


132.5 


RMS 


64.87 


Right 


337.7 


Slope 


0.3272 


Rootl 


20.29 


Root2 


243.4 





100 200 300 400 500 600 700 800 900 1000 
Mass (GeV) 



100 200 300 400 500 600 700 800 900 1000 
Mass (GeV) 



100 200 300 400 500 600 700 800 900 1000 
Mass (GeV) 



Figure 4: Invariant mass distributions for Set 1 at the SU3 point: qll thresh- 
old (left), ql low edge (middle) and ql high edge (right) 

Figures (3) and (4) show five invariant mass distributions for Set 1 of 80k 
events at the SU3 point after applied the kinematic cuts and OSOF subtrac- 
tion. For Set 1 the endpoints and reconstruction errors found from fitting of 
the edges are {103±2, 523±6, 265±4, 338±6, 461±5} to be compared with 
heoretical endpoints {103.1, 535.2, 263.8, 340.7, 456.0} for the SU3 point. 

Mass region search 

Table (4) gives the kinematic endpoints in dependence on SUSY particle 
masses [8], [9]. The different cases listed for qll and ql distributions are 
distinguished by mass ratios of sparticles. 

While the 11 edge and qll thresholds are given by unique expressions, the 
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Table 4: The kinematic endpoints for the process (3). The following notations 
are used for masses x = m |o ? I = m f > C = m |o j 9 — m | an d endpoints 

( m ZLf = (q-i)(i-t)k (KO 2 = ti-htf-x)/l (KCJ 2 = 

(q ~ — X)/(2l — x)- F° r an( i ql edges endpoints are given by different 
expressions in dependence on mass ratios. 



other endpoints have different formulas for different mass ratios which are 
not known a priori: there are four expressions for the qll edge and three 
expressions for ql low and ql high edges. Each overall combination corre- 
sponds to a unique mass region in SUSY particle mass space m g ~ L , m^o, mj- , 
m^o. These mass regions defined by mass ratios from Table (4) are labeled 
by R(i,j), where 1=1,2,3,4 and j=l,2,3 denote the corresponding mass ratios 
for qll and ql endpoints, respectively. According to the analysis performed 
in [9] only nine of these 12 (4 x 3) combinations are physical, the regions 
R(2,l), R(2,2), and R(3,3) are not possible. Additionally, the regions R(2,3), 
R(3,l), and R(3,2) are degenerate: the invariant mass distributions are not 
linearly independent, because the qll edge is the function of the 11 and ql 
high edges. 

A procedure was developed to determine from the data sample the SUSY 
particle mass-space region. Once the SUSY mass region is defined, the end- 
point expression is determined in a unique way from Table 4. Expressions for 
endpoints are inverted to get masses in terms of endpoints. Four endpoints 
are required to invert the formulas because there are four unknown masses 
in the process (3). For formula inversion we always use the 11 distribution 
because it gives the endpoint with relatively high precision and we take any 
three endpoints of the remaining four. For example, one can extract four 
sparticle masses by considering either 11, qll, ql low, ql high combination of 



endpoints or 11, qll, ql low, qll threshold combination and so on. Thus for 
each region R(l,l), R(l,2), R(l,3), R(4,l), R(4,2), R(4,3) we get four dif- 
ferent sets of mass formulas. For degenerate regions R(2,3), R(3,l), R(3,2) 
we get only two sets of mass formulas because there are only four indepen- 
dent endpoints. Note that inverse formulas are very sensitive to the endpoint 
determination precision because some formulas contain singularities in the 
range of 30 GeV around actual endpoints. The use of the full shapes of 
the signal as it was proposed in [18], [19] could lead to an improvement in 
endpoint extraction. 

The SUSY mass region is accepted if for all sets of inverse formulas the 
hierarchy condition m q ~ L > m^o > m~ lR > m^o > 20 GeV is satisfied 
and found masses correspond to the mass ratios of Table (4). Variations of 
low mass bound on the lightest neutralino mass in the range 20-50 GeV does 
not affect the results. Note that the limit from accelerator experiments on 
m^o is 46 GeV [20]. 

For all five SU3 data sets the R(l,3) region in mass space was found, 
which is the true region for the SU3 point. Note that the mass region search 
is sensitive to fluctuations in invariant mass distributions. An increase in in- 
tegrated luminosity should lead to more precise reconstruction of light SUSY 
particle masses. 

Preliminary estimation of light SUSY particle masses 

Once endpoints are found and mass regions are defined, sparticle masses can 
be determined by a minimization of x 2 ■ Parameters m of a x 2 fit are the 
four sparticle masses m q ~ L , m^o, m-j- , m^o. The x 2 ■, as a function of the free 
parameters to, is then given by 

xM)= ± (Qt^mi (4) 

i=i 17 1 

Expression Qi(m) is the endpoint as a function of sparticle masses for the cor- 
responding SUSY mass region as given in Table (4). Q° bs are the observable 
endpoints and their errors Oi determined in the above procedure. 

The results of the fit for preliminary mass estimates and their errors 
based on found endpoints, edge reconstruction errors and mass regions are 
summarized in Table (5). 

IV. Event filter 

An event filter is required before the final fit to determine the heavy SUSY 
particle (g, b) mass range and suppress background. It is very important to 
reduce combinatorial background because at the final stage of the reconstruc- 
tion procedure all possible five event combinations are considered. The com- 
binatorial background is defined as five-event combinations that include at 
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Table 5: Light masses preliminary determined from fitting of edges for the 
SU3 point. 



least one background event. The number of five-event combinations including 
only signal events is given by N s = n s !/(5!(n s — 5)!) where n s is the number of 
signal events in a data sample. The total number of five-event combinations 
including backgrounds is given by N to t = {n s + nfc)!/(5!(n s + n& — 5)!) where 
rib is the number of background events in a data sample. The contribution 
of combinatorial background with respect to all possible five signal events 
combinations (N tot — N s )/N s in the case of high n s can be approximated by 
((n s + rib)/n s ) 5 which is of about 240 for the SU3 point, as follows from Table 
3. At the event filter stage we can reduce the 5-dimensional mass space for 
each event to a 2-dimensional one, supposing that three light SUSY particle 
masses Ir, Xi) are fixed and taken as the preliminary masses found by 
the endpoint method in the previous chapter. 

The event filter procedure is based on an approximate likelihood function 
with the mass relation constraint for each event 

4 

L(mg,m~ b ) = \l exp 
i=i 

where index i runs over observed particles (two b-quarks and two leptons) 
and labels the measured absolute momenta p™ ea ^ the uncertainties <Tj in their 
measurement, and the event true absolute momenta p^ vent . The approximate 
likelihood function (5) takes into account only uncertainties in lepton and 
b-jet energy measurements. Note that positions of each of two b-jets and of 
each of two leptons in the decay chain (1) are unknown. It is quite simple 
to resolve the b-jets assignment because usually (about 96% of the time) the 
b-jet with higher pr originates from 6-quark decay, since for SU3 point at 
rest system p Tbl = (m~ — m|o)/2m^ 3> prb 2 — ( m l ~ m l)/^ m g- Therefore for 

each event we assume that the b-jet with higher p T originates from 6-quark 
decay. Leptons with higher p? are produced in both vertices with comparable 
probabilities. In this work for simplicity it is assumed that leptons with 
higher px originate from Ir decay. We use the following parametrization for 
Oi in equation (5): for b-jets a/E = 0.5/ \J E{GeV) © 0.03 and leptons 
a IE = 0.12/ y/E(GeV) © 0.005. 

To find a maximum of the likelihood function (5) with the constraint (2) 
is the same as to search for a minimum of x 2 — ~ 2 • logL function: 
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In Eq.(6) the mass relation constraint f(m,p) (2) is taking into account by 
the Lagrange multiplier A. We construct an event likelihood distribution cal- 
culating minimum of the \ 2 function (6) for each of 10 5 randomly defined 
points in the (g, b) mass plane, assuming that they are distributed uniformly 
in the range 0.4 — lATeV and 0.3 — 1.3TeV, respectively. This effectively 
corresponds to a two dimensional grid with spacing of 3GeV x 3GeV. We 
add points (nig, mi) with a weight equal to (Xcut — X 2 )/2 to an event likeli- 
hood distribution histogram if the minimization procedure of the function (6) 
converges, x 2 ( m gi m b) < Xmt = 10 and / < lO^^GeV 2 . The minimization 
iteration procedure is considered to have converged if the relative difference 
in xlvent f° r two consecutive iteration steps is less than 0.5% and the total 
number of iteration steps does not exceed 20. For signal events the event 
likelihood distribution has a maximum in the region of the (g, b) mass plane 
correlated with the true masses of g and b. Thus signal events should give a 
peak in the region of true masses. For background events there is no strong 
correlation of maximum likelihood distribution with true (g, b) masses. We 
define the combined event likelihood function as 



log L comb (mg,m 6 ) = lo £ L ( m i> m b) 



(7) 



events 



where the sum includes all events in the data sample. 

Figure (5) shows two projection histograms of combined event likelihood 
distribution versus gluino mass m g and versus difference of gluino and sbot- 



tom masses (m- n 
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for Set 1. The variables 



[m g , m- g 



mi) are used 

because an event likelihood distribution is more sensitive to these variables. 
Preliminary estimates for heavy sparticle masses and mass errors (the mean 
values and standard deviations) obtained by fitting of these histograms by a 
Gaussian near the peak region are given in Table 6 for all five sets. 
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Figure 5: Combined event likelihood projection histograms for Set 1 at the 
SU3 point. 
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649±110 


b 


569±174 


615±135 


577±162 


621±214 


536±116 



Table 6: Preliminary heavy SUSY particle masses and mass errors estimated 
by event filter. 



To suppress the background the previous procedure constructing an event 
likelihood distribution is repeated, but in the narrow range in m g and m^: 
rrig ± a, mi±a from Table 6, assuming the uniform distributions. Requiring 
convergence and applying \ 2 an d the constraint f cuts as above results in a 
likelihood distribution histogram that has high efficiency for a signal event 
and low efficiency for background events. Typically for a signal event the 
likelihood distribution histogram has of a few thousand entries (recall that 
the procedure is repeated 10 5 times). In order to reject background and 
not suppress signal, an entry threshold of 300 is chosen. Thus if the event 
likelihood distribution histogram has greater than 300 entries the event is 
considered a signal candidate and it is retained for the consequent analysis. 

After the application of the event filter the ratio of background events 
to signal events is reduced approximately by a factor 2 as can be seen from 
Table 7 . 



Set 


Total 


Signal 


SUSY 


Ratio 




Number 


Events 


Backg. 




1 


154/83 


47/39 


107/44 


3.3/2.1 


2 


131/83 


48/42 


83/41 


2.7/2.0 


3 


148/76 


44/35 


104/41 


3.4/2.2 


4 


157/110 


50/45 


107/65 


3.1/2.4 


5 


141/89 


55/45 


86/44 


2.6/2.0 


1-5 


731/441 


244/206 


487/235 


3.0/2.1 



Table 7: The number of signal and background events before/after an appli- 
cation of event filter to 80k SUSY events. The last row shows the sum over 
all five sets. 

Contribution of the combinatorial background which is given approx- 
imately by ((Signal+Background)/Signal) 5 is therefore significantly sup- 
pressed. The suppression factor for the SU3 point varies from 3 to 9. 



V. Combinatorial method for final mass recon- 
struction 

A combinatorial procedure is used for the final SUSY particle mass recon- 
struction. It is applied only to the events that pass the event filter. In order to 
explain this point let's neglect for a moment uncertainties in detected particle 
momenta and the widths of the Breit-Wigner distribution. For a single event 
the mass relation constraint represents a four dimensional surface in five- 
mass parameter space. Therefore, for two events the intersection of two four 
dimensional surfaces is a three dimensional surface and correspondingly for 
three events one gets a two dimensional surface of intersection, for four events 
the surface is one dimensional and at last for five events the intersection is 
just a point or a few points in five dimensional mass space corresponding to 
SUSY particle masses. Thus in an ideal case five events would be enough [12] 
to reconstruct masses of SUSY particles. Uncertainties in detected particle 
momenta and Breit-Wigner distribution lead to smearing in the position of 
this point. There is a trade-off between the number of events in a combi- 
nation and combinatorial background: more events restrict the number of 
degrees of freedom and better define SUSY masses, but combinatorial back- 
ground increases. Therefore, at the final stage of mass reconstruction when 
the physical background has already been reduced, we will consider all possi- 
ble five event combinations from the event sample. Recall that in reality due 
to the Breit-Wigner distribution the masses of the gluino and bottom squark 
vary from event to event. A method described below takes into account that 
SUSY particle masses can vary from event to event. 

SUSY particle masses are reconstructed by a search for a maximum of 
a combined likelihood function constructed for each possible combination of 
five events in the data sample. For sparticle masses (rh) the combined likeli- 
hood function for the combination is defined as the product of the maximum 
likelihood functions for individual events. To find a maximum of the com- 
bined likelihood function for the combination is the same as to search for a 
minimum of the function 

5 

xLj™) = E MIN (X 2 cvcnt)i (8) 
i=l 

where xLmb = ~ 2 ' logL comb . In Eq.(8) MIN(Xe ven t)i is a result of searching 
for a minimum of the xtvent function (9) for an individual event for given rh 
with the mass relation and 11 edge constraints. For each of the five events in 
the five event combination the MIN(Xe vent ) * s fitted with 9 parameters (four 
particle momenta and five SUSY masses) starting with p event = p meas and 
m event = rh. The minimization iteration procedure is converged if a relative 
difference in Xevent f° r two consecutive iteration steps is less than 0.5% and 
the total number of iteration steps does not exceed 20. 
The x 2 function for an event is defined by 



J (r^event r .meas\2 jj ( ^event „_ \2 

^ = E ~ / + E ^ } +Ai/ + A 2 /». (9) 

i=l 1 n=l ™ 

where the first term takes into account deviations of measured momenta of b- 
jets and leptons from the true ones. The second term takes into account that 
sparticles are varied from event to event and approximated by a Gaussian 
of width a n instead a Breit-Wigner distribution. In Eq.(9) the mass relation 
and 11 edge constraints are taking into account by Lagrange multiplier Ai, A 2 . 
Standard deviations corresponding to the mass term are taking to be 15 GeV 
for the gluino, 5 GeV for bottom squark and 1 GeV for light masses. The 
first two numbers are comparable with theoretical widths for heavy SUSY 
particles. The last number takes into account the fact that light SUSY 
particles are quite narrow or stable. We note that the results of the mass 
reconstruction are not strongly sensitive to actual values of sparticle widths. 

To get a starting point for the minimization of Xcomb f° r eac h five-event 
combination we calculate Xcomb f° r 3000 m values produced with a simple 
Monte Carlo. The set of the mass rh with the smallest Xcomb * s used as the 
starting point. It is assumed that heavy masses are distributed uniformly in 
the following range: mean value ± 2a, where the mean values and standard 
deviations are given in Table 6 as the results of event filter application. 
For the three light masses the Gaussian distribution is assumed with mean 
values and standard deviations found by the endpoint technique as given 
in Table 5. The starting point is used by MINUIT code with the Simplex 
algorithm to search for the minimum of Xcomb w hh five parameters rh. The 
five-event combination is added to reconstruction mass histograms if the 
MINUIT minimization procedure converges, Xcomb < Xcut = 10 and the sum 
of the mass relation and 11 edge constraint functions for five events is less than 
1 GeV 2 . The CPU time required for the final minimization of the five-event 
combination is about 2 second for a 3GHz processor. For a set of 100 events 
the required time for combinatorial mass reconstruction would be about 7 
years. In order to reduce the computational time, we divide a sample of 
80k SUSY events into four subsets of 20k. Forming combinations only from 
events within each 20k subset greatly reduces combinatorics. Because the 
total set is subdivided in four subsets the procedure is carried on for each 
subset and reconstruction mass histograms are merged. 

Reconstructed SUSY particle mass distributions and results of fitting 
these distributions by a Gaussian are presented in Figures (6 and 7) for 
Set 1 at the SU3 point with integrated luminosity 4.2 fb _1 . As can be seen 
in these Figures, the reconstructed mass distributions are well described by 
a Gaussian. The mean values of the Gaussian fit are considered as the re- 
constructed sparticle masses. As can been seen these masses are close to the 
theoretical masses. It is assumed that mass reconstruction errors are given 
by the Gaussian fit standard deviations. For the light sparticles the final pro- 



cedure improves mass reconstruction errors by a factor of 2 but the masses 
are changed only slightly in comparison with the preliminary estimate. 




Figure 6: Reconstructed heavy SUSY particle mass distributions for Set 1, 
including background with integrated luminosity 4.2 fb _1 : gluino g (left), 
sbottom b (right). The line is a result of Gaussian fit. 




Figure 7: Reconstructed light SUSY particle mass distributions for Set 1, in- 
cluding background with integrated luminosity 4.2 fb -1 : neutralino xX (left); 
slepton l R (middle), neutralino Xi (right). The line is a result of Gaussian 
fit. 

Final results of this mass reconstruction approach are presented for five 
data sample sets of 80k events each in Table 8. 

In order to illustrate a spread in reconstructed masses the results of Table 
8 are also shown in a form of ideograms [21] in Figures (8 and 9) for the five 
data sample sets. Each reconstructed mass in an ideogram is represented by 



Set 


m~ g 


m- h 


m x° 


m l 




1 


701±57 


600±57 


208±21 


148±23 


98±22 


2 


712±55 


608±53 


254±21 


183±23 


143±20 


3 


664±78 


564±80 


219±24 


149±24 


109±23 


4 


767±62 


649±65 


258±35 


193±35 


148±34 


5 


655±45 


545±47 


208±21 


138±22 


96±20 



Table 8: Reconstructed SUSY particle masses and reconstruction errors for 
five data sample sets of 80k events each. 



a Gaussian with a central value rrii, error a, and area proportional to l/<7j 
The solid curve is a sum of these Gaussians. 
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Figure 8: Ideograms of reconstructed heavy SUSY particle masses, including 
background for five data sample sets with integrated luminosity 4.2 fb : 
gluino g (left), sbottom b (right). The triangle marker gives the position of 
theoretical mass. Points with error bars correspond to five data sample sets. 



The Gaussian form of ideograms and relatively small shift of peak posi- 
tions with respect to theoretical masses demonstrate the self-consistency of 
the mass reconstruction approach. 



IX. Conclusion 

We have developed an approach that allows extracting all SUSY particle 
masses for long cascades at the LHC with an integrated luminosity of a 
few fb -1 . This luminosity can be reached at the early stage of the LHC in 
comparison with a projected integrated luminosity of 300 fb -1 in five years. 
This approach is based on a consecutive use of the endpoint method, an event 
filter and a combinatorial mass reconstruction method. 



| Neutralino /" Mass | 






SleptonfMass | 






| Neutralino x° Mass 
































HI II 1 1 II II 1 1 1 II II II 


III! 


II II II 






nil 


I i , 


III III i n III i ll 










Mi III II 


III 




















































▼■■ 






































































































































































































































III 


i i i i 


I llU 


', II , 


I I I I 


, , , , i\ 


i i i i 


, J, \ , , . , i , 








i , J 


i i i i 


i i i i 


i i i i 




1 1 1 1 







50 100 150 200 250 300 
Mass (GeV) 




c 


50 100 150 200 250 300 
Mass (GeV) 




C 


50 100 150 200 250 300 
Mass (GeV) 





Figure 9: Ideogram of reconstructed light SUSY particle mass distributions, 
including background for five data sample sets with integrated luminosity 
4.2 fb _1 : neutralino X2 (left), slepton l R (middle), neutralino Xi (right). 
Triangle marker gives the position of theoretical mass. Points with error 
bars correspond to five data sample sets. 

The endpoint method allows preliminary estimate of light SUSY particle 
Ir-> Xi) masses in a model-independent way. 

The light sparticle masses are used to construct the maximum likelihood 
distribution in the two heaviest-sparticle mass-plane for each individual event 
taking into account the mass relation equation constraint. The event filter 
combining the maximum likelihood distributions for all events in the data 
sample allows a determination of the range of heavy sparticle masses and sig- 
nificant suppression of background. This is a new technique for background 
suppression. Note that in this work we did not use r tagging for background 
suppression. The application of r tagging could provide an additional sup- 
pression of combinatorial SUSY background by approximately a factor of 
2. 

Finally, SUSY particle (g, b, x°, Ir, x?) masses are extracted by a search for 
a maximum of a combined likelihood function constructed for each possible 
combination of five events in the data sample satisfying the filter. Mass peaks 
for five sparticles are clearly reconstructed because the remaining background 
does not exhibit peaks corresponding to the signal region. With more events 
the same technique could be used to reconstruct, for example, masses of the 
fx and ti in the following cascade decays: g — > b — > X2 ~^ ^1 > Xi an d 

9 -»• h -»• xl -»• Ir -» Xi- 

SUSY data sample sets for this work were generated for the benchmark 
SU3 point of mSUGRA scenario. A detector response was parametrized by 
the AcerDET detector simulator. The stability and precision of the approach 
was demonstrated by considering five different 80k event data sets for each 
model point. Masses were reconstructed with a precision of about 10% for 



heavy sparticles and 10%, 15% and 20% for light yQi an d Xi sparticles, 
respectively. The precision of mass reconstruction with this technique should 
be improved with an increase in integrated luminosity. 

We also applied the cascade mass reconstruction approach to SPSla 
model point of mSUGRA parameter space and masses were reconstructed 
with the same precision as for the SU3 point. 

The approach developed in this paper can be used for the ATLAS and 
CMS detectors at the LHC; or at the future ILC. 
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Appendix A. Mass relation equation 



The mass relation method [12] is based on a solution of a system of kinematic 
equations obtained for each vertex of decay cascade. For the chain (1) the 
system of kinematic equations is given by four equations corresponding to 
four vertices 

\2 



m 2 ~ 
x 2 

m 



b 



(Pxi + k h + h 2 ) 2 

(p x ° + k h + k h +PbJ 2 
(P*° + k h + k h + Pbi +Pb 2 f 



This system of four equations can be easily modified to a linear system 
defining four component of lightest neutralino four-momentum 



where matrix S is given by 



Sp x o = Q, 



(10) 



( k h k h k h k h\ 

h° —h l —h 2 —k- 3 

K h K h K h K h 

P° bl -Pi ~P\ -P 3 bl 



\p° b2 -p1 2 -p\ 2 -p\J 



p^o is the column vector of lightest neutralino four-momentum 



Q is the column vector of coefficients 



where 



c-k h - k h 
d-Pbi- (kh + h 2 ) 
\e-p b2 ■ (k h + k h +p hl )J 



b = 0.5 * (m? — m 2 ~ ) 
c = 0.5 * (m? — m- ) 

V X 2 \r> 

m|o) 



d = 0.5 * (ml 

v b 



0.5 * (m? 



Lepton and quark masses are neglected in these equations. 

The solutions of linear equation (10) is given by standard formulas 



pl a = detSj/detS 



(11) 



where the submatrix Sj is formed by substituting elements of vector Q 
instead of j th column of matrix S. 

By using on-shell mass condition m~ = p 2 ~ for the lightest neutralino 

Xi Xi 

with momentum components from (11) one gets the following equation 

j=3 



Equation (12) is the basic equation of mass relation method and it gives 
the relation between masses of five SUSY particles. Each event is represented 
as a curve in five dimensional mass space. Coefficients of this equation are 
functions of the four-momenta of detected particles where b quarks are mea- 
sured as jets in the detector. 




(12) 



